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The baryon distribution is studied by using a parton cascade model which is based on pQCD 
incorporating hard partonic scattering and dynamical hadronization scheme. In order to study 
baryon distribution, baryonic cluster formation is newly implemented as well as hadronic higher 
resonance states from parton/beam cluster decay. The net baryon number and charged hadron 
distributions are calculated with different TVfactors in which parameters are fixed by elementary pp 
data at E c . m . = 200 GeV. It is found that baryon stopping behavior at SPS and RHIC energies are 
not consequence of hard parton scattering but soft processes. 
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I. INTRODUCTION 

Heavy ion experiments at BNL-AGS and CERN-SPS have been performed motivating by the possible creation 
of QCD phase transition and vast body of systematic data such as proton, pion strangeness particles distributions, 
HBT correlation, flow, dileptons and J /if) distributions have been accumulated including mass dependence and their 
excitation functions Jl|-^| ■ Data from forthcoming experiment at BNL- RHIC will be available soon. 

Strong stopping of nuclei has been reported both at AGS and at SPS energies It is reported that baryon 

stopping power can be understood within a hadronic models if we consider multiple scattering of nucleon using 
reasonable pp energy loss ||. For example, within string based models |7[-pT|, baryon stopping behavior at SPS 
energies is well explained by introducing diquark breaking mechanism in which diquark sitting at the end of the string 
breaks. Diquark breaking leads to large rapidity shifts of the baryon. Constituent quark scattering within a formation 
time H[l2| has to be considered in order to generate Glauber type multiple collision at initial stage of nuclear collisions 
in microscopic transport models which describe full space-time evolution of particles. 

Event generators based on perturbative QCD (pQCD) are proposed such as HIJING (Heavy Ion Jet Interaction 
Generator) (l^Jl^], VNI (Vincent Le Cucurullo Con Giginello) Jl5|], in order to describe ultra-relativistic heavy ion 
collisions emphasizing the importance of mini-jet productions. VNI can follow the space-time history of partons and 
hadrons. The parton cascade model of VNI has been applied to study several aspects of heavy-ion collisions even 
at SPS energies [p^| . However, original version of VNI implicitly assumed the baryon free region at mid-rapidity 
during the formation of hadrons, because only two parton cluster (mesonic cluster) formations are included in the 
Monte-Carlo event generator VNI 

In this work, The baryon distribution at SPS and RHIC energy are discussed using modified version of parton cascade 
simulation code VNI fl6|| . The main features of the parton cascade model to be used here are that implementation 
of baryonic cluster formation and during the parton/beam cluster decay higher hadronic resonance states are allowed 
to produce in order to be able to calculate baryon distribution in heavy ion collisions. 



II. PARTON CASCADE MODEL 



First of all, the main features of the parton cascade model of VNI as well as the main points of the modification 
will be presented. Relativistic transport equations for partons based on QCD [ fi8| are basic equations which are 
solved on the computer in parton cascade model. The hadronization mechanism is described in terms of dynamical 
parton- hadron conversion model of Ellis and Geiger p9|-pl||. The main features in the Monte Carlo procedure are 
summarized as follows. 

1) The initial longitudinal momenta of the partons are sampled according to the measured nucleon structure function 
f(x, Qq) with initial resolution scale Qq. We take GRV94LO (Lowest order fit) |^2| for the nucleon structure function. 
The primordial transverse momenta of partons are generated according to the Gaussian distribution with mean value 
of pj_ = 0.44GeV. The individual nucleons are assigned positions according to a Fermi distribution for nuclei and the 
positions of partons are distributed around the centers of their mother nucleons with an exponential distribution with 
a mean square radius of 0.81fm. 
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2) With the above construction of the initial state, the parton cascading development proceeds. Parton scattering 
are simulated using closest distance approach method in which parton-parton two-body collision will take place if their 
impact parameter becomes less than y^a/ir, where a represents the parton-parton scattering cross section calculated 
by pQCD within a Born approximation. Both spacelike and timelike radiation corrections are included within the 
leading logarithmic approximation. Elementary 2-^2 scatterings, 1—^2 emissions and 2 — > 1 fusions are included in 
the parton cascading. 

3) Parton clusters are formed from secondary partons that have been produced by the hard interaction and parton 
branching. The probability of the parton coalescence to form color-neutral cluster n is defined as ]2(J 

0, Lij < Lq, 



Lq — 



Uij^c = { 1 - exp( L _ L -X Lo < Uj < L c , (1) 
lj Lij > L c , 

where L c = 0.8fm is the value for the confinement length scale and Lq = 0.6fm is introduced to account for finite 
transition region. L,j is defined by the distance between parton i and its nearest neighbor j: 

L i:j ee min(Aji, • • • , Ay, • • • , A m ), (2) 



where Ay = */ rfrj^ is the Lorenz-invariant distance between partons. So far, only the following two-parton coales- 
cence 

g + g^Ci + C 2 ,g + g^C + g,g + g^C + g + g, (3) 
q + q^C\+C 2 ,q + q^C + g, (4) 
q + g^C + q,q + g->C + g + q. (5) 

have been considered in the VNI model. In this work, if diquarks are formed with the above formation probability, 
baryonic cluster formation is included as 

qq + q^C, (6) 
qq + q^C, (7) 
9l92 + 93 -> 9l93 + 92, (8) 
9l92 +g — ► 919293 + 93- (9) 

Note that by introducing those cluster formation processes, We do not introduce any new parameters into the model. 

4) Beam clusters are formed from primary partons (remnant partons) which do not interact during the evolution 
even though they travel in the overlapping region of nuclei. They may be considered as the coherent relics of the 
original hadron wavefunctions, and should have had soft interactions. Those underlying soft interactions are simulated 
by the beam cluster decay into hadrons in VNI because additional possibility that several parton pairs undergo soft 
interactions. This may give a non-negligible contribution to the 'underlying event structure' even at the collider 
energies. The primary partons are grouped together to form a massive beam cluster with its four-momentum given 
by the sum of the parton momenta and its position given by the 3- vector mean of the partons' positions. 

5) The decay probability density of each parton cluster into final state hadrons including hadronic resonances is 
chosen to be a Hagedorn density state. The appropriate spin, flavor, and phase-space factors are also taken into 
account. In the decay of parton/beam cluster, higher hadronic resonance states up to mass of 2GeV can be produced 
in our model. 

To summarize, the main different points from original version are 1) baryonic cluster formation. 2) inclusion of 
higher hadronic resonance up to mass of 2GeV. 3) exact conservation of flavor, i.e. (baryon number, charge, etc). 
4) reasonable total momentum conservation: total momentum is conserved within 10% at RHIC energy for central 
Au+Au collision. 

III. RESULTS 
A. Elementary collisions 

Since our version of parton cascade code differs from original version of VNI, we have to check the model parameters. 
First, particle spectra from pp collisions at y/s = 200GeV calculated by the modified version of VNI are studied to 
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see the model parameter dependence. Here we see the K- factor dependence as mentioned in Ref. [ p3| . In Fig. [|, 
experimental data on pseudorapidity distributions (left panel) and the invariant cross sections (right panel) are 
compared to the calculation of the parton cascade model with different parameters on the treatment of so-called 
if-factor. The calculations (upper three figures) are done by adding the constant factor to the reading-order pQCD 
cross sections: 

<?pQCd(Q 2 )=K xa LO (Q 2 ) 

with values K — 1,2,2.5. While bottom figure corresponds to the calculation changing the Q 2 scale in the running 
coupling constant a s as 

& p qcd(Q 2 ) = v LO {a s {vQ 2 )) 

with the value r/ = 0.075. We also plot the contribution from parton cluster decay in the left panel with dotted 
lines. The contribution of parton cluster decay which is come from interacted parton coalescence changes according 
to the choice of the correction scheme. We can fit the pp data of pseudorapidity distributions with different correction 
schemes as seen in Fig. [I] by changing the parameter (in actual code, parv(91)) which controls the multiplicity from 
beam cluster. We have to check the model with various elementary data including incident energy dependence in 
order to fix model parameters. Next we will present some results on nuclear collisions with those parameters. 





charged hadrons, lr|l<2.5 
ti= 0.075 




2 4 6 

p_L (GeV/c) 



3 



FIG. 1. Data on charged particle ((h~ + h + )/2) pseudorapidity (left) and the invariant cross sections Ed 3 a/d 3 p 
(right) ]2^] from pp collision at ^fs = 200GeV compared to parton cascade model calculations with various parameters. The 
contributions from parton cluster decay are also plotted by dotted lines in a left panel. 



B. Comparison with SPS data 

The baryon stopping problem is one of the important element in nucleus- nucleus collisions. Original version of 
VNI implicitly assumed baryon free region at midrapidity, because baryonic parton cluster formation is not included. 
Baryons only come from beam cluster, not parton cluster formation in the original version of VNI. We can now discuss 
the baryon stopping problem with our modified version of VNI. 

We have calculated the net proton distribution at SPS energy to show the reliability of the modeling of beam cluster 
formation in the parton cascade model. Fig. || compares the parton cascade calculation for Pb+Pb collision at the 
laboratory energy of Ei a t, — 158 AGeV with the if-factor 1.0 (original version uses rj — 0.035) of net protons with the 
data |||. It is seen that contribution from parton cluster is neglibigly small, thus baryon stopping behavior is fully 
explained by soft physics (in this case, beam cluster decay) when we chose the if- factor 1.0 at SPS energies. It should 
be noted that there is no microscopic dynamics in the modeling of the beam cluster formation in the parton cascade 
model, but it is a simple fit to the data of pp collisions. 
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FIG. 2. Parton cascade model calculations of the rapidity distributions of net protons for Pb + Pb collision at SPS energy 
(Eiab = 158 AGeV). K — 1.0 (solid and dotted lines) and K = 2.0 (dashed and dash-dotted lines) are used in this calculation. 
Dotted and dash-dotted lines corresponds to the contribution from parton cluster decay respectively. 



C. Predictions for RHIC 



The if-factor dependence of both net proton and charged particle rapidity distribution are studied in Fig. || In 
terms of net proton distribution, there is no strong if-factor dependence. We can see that parton cluster formation 
and its decay predict almost baryon free at mid-rapidity region regardless of the choice of if-factor, though there are 
lots of protons and antiprotons at mid-rapidity. We conclude that hard parton scattering plays no rule for the baryon 
stopping within a parton cascade model. However, note that string based model like HIJING/B |) 10 predicts proton 
rapidity density of 10 and UrQMD predicts [^6| 12.5 at mid-rapidity. However, as pointed out in Ref. 0), charged 
hadron multiplicity is strongly depend on how to chose the leading order correction scheme. 
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FIG. 3. Parton cascade model calculations of the rapidity distributions of net protons, protons, and anti-protons (upper) 
and charged particles (lower) for Au + Au collision at E c , m , = 200 AGeV for head on collisions. 




Fig. ^ displays the net baryon number distributions as a function of rapidity obtained by parton distribution from 
parton cascade before hadronization with the if-factor of 1 (left) and 2.5 (right). Net baryon number of time-like 
partons are distributed around mid rapidity region but its contribution are small as consistent with the net proton 
distribution in Fig. |^. 
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FIG. 4. Rapidity distributions of net baryon number (q — q) obtained from parton cascade model before hadronization with 
K = 1 (left) and K = 2.5 (right) for Au + Au collision at E c . m . = 200 AGeV for head on collisions. 
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IV. SUMMARY 



In summary, first, we have checked that different treatments for the inclusion of higher-order pQCD corrections 
in parton cascade model can fit the elementary pp collisions. We have to check other elementary processes to fix 
the model parameters. We show the net proton rapidity distribution at SPS energies to demonstrate that the beam 
cluster treats underlying soft physics in the parton cascade model reasonably well for nucleus nucleus collisions. Then, 
we have calculated the net proton rapidity distribution at RHIC energy as well as charged particle distributions using 
modified version of parton cascade code VNI in which we newly introduced baryonic parton cluster formation and 
higher hadronic resonance states from decay of parton and beam cluster. Within a framework of perturbative parton 
cascading and dynamical hadronization scheme, we predict almost baryon free plasma at RHIC energy. The charged 
particle rapidity distributions are also studied with the parameter set which are fitted by pp collisions. Strong if- factor 
dependence on the hadron multiplicity is seen as previously being found by Ref. [ p3[ . we can not fix the if- factor 
from only rapidity and transverse momentum distributions for pp collisions. 

In this work, we consider only two or three parton coalescence, but in dense parton matter produced in heavy ion 
collisions, this assumption might be broken down. Inverse processes like hadron conversion to parton such as C — ► qq 
are also ignored which might become important at higher colliding energies. 
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